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Abstract 

Bedau has developed a general set of evolutionary statis- 
tics that quantify the adaptive component of evolution- 
ary processes. On the basis of these measures, he has 
proposed a set of 4 classes of evolutionary system. All 
artificial life sytems so far looked at fall into the first 3 
classes, whereas the biosphere, and possibly the human 
economy belongs to the 4th class. The challenge to the 
artificial life community is to identify exactly what is 
difference between these natural evolutionary systems, 
and existing artificial life systems. 

At ALife VII, I presented a study using an artificial evo- 
lutionary ecology called ^ c tab. Bedau's statistics cap- 
tured the qualitative behaviour of the model. E- c< Lab ex- 
hibited behaviour from the first 3 classes, but not class 
4, which is characterised by unbounded growth in di- 
versity. E c H.ab exhibits a critical surface given by an 
inverse relationship between connectivity and diversity, 
above which the model cannot tarry long. Thus in or- 
der to get unbounded diversity increase, there needs to 
be a corresponding connectivity reducing (or food web 
pruning) process. This paper reexamines this question 
in light of two possible processes that reduce ecosystem 
connectivity: a tendency for specialisation and increase 
in biogeographic zones through continental drift. 



Introduction 

During the Phanerozoic (540Mya-present), the diversity 
of the biosphere (total number of species, also known 
as biodiversity) has increased dramatically. The trend 
is most clear for intermediate taxonomic levels (fam- 
ilies and orders), as fossil species data is too incom- 
plete and higher taxonomic levels (phylum and class di- 
versity) have been fairly constant since the Paleozoic. 
A recent review is given by Benton(2001). The most 
completely documented diversity trend is amongst ma- 
rine animals, which exhibits a plateau during the Paleo- 
zoic (540-300Mya), followed by an accelerating diversity 
curve since the end of the Permian. The correspond- 
ing trend amongst continental, or land animals is char- 
acterised by a clear exponential growth since the first 
species colonised dry land during the Ordovician. Ben- 
ton argues that the terrestrial trend is more character- 
istic than the marine trend, owing to the far greater 



diversity shown amongst land animals, even though the 
marine fossil record is more complete. Similar trends 
have been reported for plants (Tiffney & Niklas 1990). 

Bedau et aL(1998) introduced a couple of measures to 
capture the amount of adaptation happening in a general 
evolutionary system. The basic idea is to compare the 
dynamics of the system with a neutral shadow system 
in which adaptation is destroyed by randomly mixing 
adaptive benefits amongst the components of the sys- 
tem (think of the effects ultra-Marxism might have on 
an economy!). The amount of adaptive activity (num- 
bers of each component in excess of the shadow model 
integrated over time) and adaptive creativity (numbers 
of speciations per unit time exceeding a threshold of ac- 
tivity) is measured. Bedau has also introduced a general 
neutral shadow model that obviates the need to generate 
one on a case by case basis (Rechtsteiner & Bedau 1999). 

Using these measures, it is possible to distinguish 3 
classes of activity: 

1. unadaptive evolution, when the mutation rate is so 
high that organisms have insufficient time to have their 
adaption tested before being killed off by another mu- 
tation 

2. adapted but uncreative evolution, when species are 
highly adapted, but mutation is so low that ecosystems 
remain in perpetual equilibrium 

3. creative, adaptive evolution, when new species contin- 
uously enter the system, and undergo natural selection 

The Biosphere appears to be generating open ended 
novelty — not only is it creative, but it is unbound- 
edly creative. Evidence for this exists in the form of 
the intricate variety of mechanisms with which differ- 
ent organisms interact with each other and the environ- 
ment, and also in the sheer diversity of species on the 
planet. Whilst there is no clear trend to increasing or- 
ganismal complexity (McShea 1996), there is the clear 
trend to increasing diversity mentioned above, which is 
likely to be correlated with ecosystem complexity. Be- 
dau takes diversity as a third evolutionary measure, and 
distinguishes between bounded and unbounded creative 



evolution, according to whether diversity is bounded 
or not. All artificial evolutionary systems examined to 
date have, when creative, exhibited bounded behaviour 1 . 
This was also the case of the E- c< Lab model (Standish 
2000). Bedau has laid down a challenge to the artificial 
life community to create an unbounded, creative evolu- 
tionary system. 

Ecosystem Complexity 

The heart of the idea of unbounded creative evolution- 
ary activity is the creation and storage of information. 
The natural measure of this process is information based 
complexity, which is defined in the most general form 
in(Standish 2001). The notion, drawing upon Shannon 
entropy and Kolmogorov complcxity(Li & Vitanyi 1997) 
is as follows: 

A language L\ = (S, /i), is a countable set of possible 
descriptions S, and a map {0, 1}. We say that 

s, s' € S have the same meaning iff /z(s, s') = 1. Denote 
the length of s as £(s) and S n = {s G S : £(s) = n}. The 
information content (or complexity) of a description s is 
given by: 

card({s' e S n : m(s,s') = 1}) , v , 

C {s = - hm log 2 w g ; (1) 

n->oo card(S n ) 

In the usual case where the interpreter (which defines /j.) 
only examines a finite number of symbols to determine a 
string's meaning, C(s) is bounded above by £(s) log 2 B 
where B is the size of the alphabet. This is equivalent 
to the notion of prefix codes in algorithmic information 
theory 

Now consider how one might measure the complexity 
of an ecosystem. Diversity is like a count of the number 
of parts of a system — it is similar to measuring the 
complexity of a motor car by counting the number of 
parts that make it up. But then a junkyard of car parts 
has the same complexity as the car that might be built 
from the parts. In the case of ecosystems, we expect the 
interactions between species to be essential information 
that should be recorded in the complexity measure. But 
a simple naive counting of food web connections is also 
problematic, since how do we know which connections 
are significant to a functioning ecology? 

To put the matter on a more systematic footing, con- 
sider a tolerance e such that x, y G R are considered 
identical if \x — y\ < e. Now two different population 
dynamics x — f(x) and x — f'(x), where 

x e M™+ = {x e R" : x t > 0}, 

can be considered identical (i.e. fi(f, /') = 1 iff 2 

!/(*)- /'(aOloo <e, VxeM"+. (2) 

Channon(2001) claims his Geb artificial life system ex- 
hibits unbounded creative behaviour 

2 As an anonymous referee pointed out, trajectories dec- 
sribed by / and /' may diverge exponentially in time, and 



At this point for the sake of concreteness, let us con- 
sider Lotka-Volterra dynamics: 

x = r*x + x*/3x (3) 

where * refers to elementwise multiplication, r is the net 
population growth rate and [3 is the matrix of interspe- 
cific interaction terms. 

Over evolutionary time, the growth coefficients rj, the 
self-interaction coefficients [3u and the interspecific in- 
teraction coefficients (3ij ,i ^ j form particular statistical 
distributions p r (ri),Pd{Pu) and Po(Aj) repectively. 

Since inequality (2) must hold over all of the positive 
cone R ra+ , it must hold for population density vectors 
|a;| <C 1 and \x\ » 1. In which case eq. (2) can be 
broken into independent component conditions on r and 
(3 can be written: 

Iri-r'il < e, Vz (4) 
ll/3-/3'Hoo < e. (5) 

Since these conditions are independent, they con- 
tribute additively to the overall complexity (1). The 
term for the growth coefficients is given by: 

C r = -log 2 n/ PrWdrl 

w -^log 2 Pr(r») -Dlog 2 2£ (6) 

where eCl, and D is the ecosystem diversity. 

The complexity term for the interaction terms is given 

by 

cp = -io g2 / nPoip'^npM^iid^ 

~ l0 S2 Po {Pij ) + l0 S2 Pd + 

i=jLj i 

^ 2 log 2 2e-l (7) 

If e is chosen very small, the total ecosystem complex- 
ity is proportional to D 2 . This is because the zeros of the 
interaction matrix are encoding information. However, 
if e = 2p„(o) ' tn cn (?) becomes: 

C p = D 2 C{\og 2 p o {l3 l0 )) + D(log 2Pd (p u )) + o(D) (8) 

This gives flesh to our intuitive notion that complexity 
should somehow be proportional to the number of con- 
nections making up the food web. 

Empirically, Lotka-Volterra dynamics has been shown 
to exhibit an inverse relationship between connectivity 

that a better definition of equivalence would also require sim- 
ilarity of the attractor sets as well. The results derived here 
would only be a lower bound of the ecosystem complexity 
under this more refined definition of equivalence. 



and diversity D cx C _1 (Standish 1998). May(1972) 
demonstrated this relationship in connection with dy- 
namical stability. However, it seems unlikely that an 
ecology undergoing evolution is often stable. If this 
result holds more generally, it implies that complex- 
ity is directly proportional to diversity, so that diver- 
sity indeed is a good proxy for ecosystem complex- 
ity Although earlier foodweb studies demonstrated this 
hyperbolic diversity-connectivity relationship, more re- 
cently collected data suggests a relationship of the form 
D cx C~ 1+e , with e wO.3-0.4 (Drossel & McKane 2002). 
If complexity indeeds scales superlinearly with diversity 
as suggested by latter data, then a system displaying 
open-ended diversity growth is indeed growing in com- 
plexity, however a system displaying bounded diversity 
growth may still be growing in complexity. 

Ec tab 

is an evolutionary ecology, and is to my knowl- 
edge the first published account of population dynam- 
ics being linked to an evolutionary algorithm(Standish 
1994). The next model to be developed in this genre 
is Webworld(Drossel, Higgs, & McKane 2001), which 
features a more realistic ecological dynamics, and han- 
dles resource flow issues better. Other models in this 
genre have appeared recently (Christensen et al. 2002; 
Anastasoff 2000). 

E c 1Lab is also the name of a software pack- 
age used for implementing this model, as well as 
other models. The software is available from 
http://parallcl.hpc.unsw.edu.au/rks/ecolab. 

The model consists of Lotka-Volterra ecology: 

h = r*n + n* (3n + mutate(/x, r, n) + 7 * V 2 n. (9) 

n is the population density vector, r the growth rates 
(net births-deaths in absence of competition), (3 the in- 
teraction matrix, n the (species specific) mutation rates 
and 7 the migration rate. In the panmictic case, the 7 
term is left out, and n refers to total populations, rather 
than population densities. 

The mutation operator randomly adds new species i 
into the system with phenotypic parameters {r il ftij, fii 
and 7j) varied randomly from their parent species. A 
precise documentation of the mutation operator can be 
found in the ^ c £a£> Technical Report(Standish ). 

The n vector has integral valued components — in as- 
signing a real valued vector to it, the values are rounded 
up randomly with probability equal to the fractional 
part. For instance, the value 2.3 has a 30% probabil- 
ity of being rounded up to 3, and a 70% probability of 
being rounded down. Negative values are converted to 
zero. If a species population falls to zero, it is considered 
extinct, and is removed from the system. It should be 
pointed out that this is a distinctly different mechanism 



than the threshold method usually employed to deter- 
mine extinction, but is believed to be largely equivalent. 

Diversity D is then simply the number of species with 
rii > 0, and connectivity is the proportion of interspecific 
connections out of all possible connections: 



Spatial Ec<n a b [ s implemented as a spatial grid, with 
the V 2 term being replaced by the usual 5-point stencil. 

Specialisation 

A specialist is a species that only depends on a restricted 
range of food sources, as opposed to a generalist which 
might depend on many food sources. A specialist has 
fewer incoming predator-prey links in the food web than 
does a generalist. Much evolutionary variety is expressed 
in sophisticated defence mechanisms that serve to sup- 
press outgoing predator-prey links. In this context, I will 
use the term specialist in a more general sense to refer 
to species with a small number of food web links. In 
order for the panmictic EcoLab model to generate an in- 
creasing diversity trend, a corresponding specialisation 
trend must also be present (which it isn't in the case 
of the usual mutation operator). Interestingly, special- 
isation is usually considered to be the default mode of 
evolution(Vermeij 1987). Generalists only exist because 
they happen to be more robust against environmental 
perturbation. 

This experiment involves modifying the mutation op- 
erator to bias it towards removing interaction terms. 
The usual ^ Cl £a£> mutation operator operator adds or 
removes connections according to |l/ r Ji where r G 
(— 1, 1) is a uniform random variate(^ c £ab Technical 
Report). In this experiment, a new experimental pa- 
rameter g G (—1, 1) (gen_bias) is introduced such that 
r G (— 1 +g, 1 + g) and the number of connections added 
or deleted is given by [(1 + s gn(r)g)/rJ ■ By specifying a 
very negative value of g, the mutation operator will tend 
to produce specialists more often than generalists. The 
code for this experiment is released as ^ c £ab 4.2. 

A typical run with g = —0.9 is shown in Figs. 1-4. 
As described in (Standish 2000), activity is weighted by 
the population density, not just presence of a particu- 
lar species. The results show unbounded creative evo- 
lutionary activity (Class 4 behaviour). As can be seen 
from Fig. 4, the system remains close to the hyperbolic 
critical surface, yet the dynamic balance has been re- 
moved by the specialisation trend. If we assume that 
D < Co/C(D), then unbounded diversity growth can 
only happen if C vanishes at least as fast as 1/D (see 
Fig. 5). An ecosystem consisting entirely of specialists 
has a constant number of foodweb links per species, or 
C cx 1/D. The presence of generalists in the ecosystem 
damps the growth in diversity, and unbounded growth is 
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Figure 1: Diversity growth for a typical run with g = 
-0.9 





Figure 4: Connectivity vs Diversity for the same run 
depicted in Figs. 1-3. The dashed curve is 8/D. 
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Figure 2: Mean cumulative activity A cum (t) 
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Figure 5: Diversity is constrained to lie under the curve 
1/C(D). The intersection of this curve with line y = D 
gives the maximum possible diversity in the ecosystem. 
If C(D) < o(D~ 1 ), then diversity is unbounded. 
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Figure 3: Cumulative New Activity over time 



only possible if the proportion of gcneralists continually 
diminishes over time. 

Continental Drift 

In (Standish 2000), I suggested that one possible expla- 
nation for the diversity growth since the end of the Per- 
mian was the breakup of the supercontinent Pangaea. 
A simple estimate given in that paper indicated that 
the effect might account for a diversity growth of about 
3.5 times that existing during the Permian. This was 
remarkably similar to the growth reported by(Benton 
1995), however it is worth noting that Benton's data re- 
ferred to families, not species. It is expected that the 
numbers of species per family also increased during that 
timc(Benton 2001). Furthermore, when continental or- 
ganism are included, familial diversity today is more like 
5 times the diversity during the Permian. 

Unbeknownst to me at the time, Vallcntinc(1973) 
had proposed essentially the same theory, called biogeo- 
graphic provincialism (the notion that the number of bi- 
ological provinces is increased through rearrangement of 
the continents). The idea received some serious support 
by Signor(1990), although in a later review he was less 
enthusiastic(Signor 1994). Tiffney and Niklas(1990) ex- 
amined plant diversity in the northern hemisphere and 
concluded that plant diversity correlated more with the 
land area of lowlands and uplands, rather than continen- 
tal breakup. Benton(1990) is characteristically sceptical 
of biogeographic provicialism as an explanation of the 
diversity trend through the Phanerozoic. Biogeography 
theory depends on an assumed dynamic balance between 
speciation and extinction 3 , which appears to be contra- 
dicted by the fossil data for continental animals(Benton 
2001), which shows a strong exponential increase in di- 
versity through the Phanerozoic. 

Since the ^ c tab model has this dynamic balance be- 
tween speciation and extinction when the dynamics self- 
organise to the critical surface D « (Cs 2 ) -1 , I exper- 
imented with the spatial version of ^ c tab reported in 
(Standish 1998). The maximum migration rate I7I00 
was swept up and down exponentially in time accord- 
ing to 0.9*/ 1000 , i.e. with a time constant of about 
9500 timesteps, by scaling 7 by 0.9 every 1000 timesteps 
(and then inverting the scaling factor every 174,000 
timesteps). It is a little hard to relate ^ c tab figures 
to biological evolution. The maximum growth rate in 
Ec<t a b is 0.01, so the doubling time for the fastest or- 
ganism in the ecosystem is around 100 timesteps. This 
might correspond to a year or so of real time. So migra- 
tion rates are being forced much faster than is typical 
in the real world. However, in EcoLab we also tend run 

3 Benton calls this a dynamic equilibrium, although it is 
nothing like what the term equilibrium denotes in dynami- 
cal system theory, and characterises biogeograhic theories as 
equilibrial 
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Figure 6: Typical run sweeping the maximum migration 
rate |7|oo- The migration operator was applied every 
100 time steps, so the units of the x-axis are 0.01 cells 
per timestep 

the mutation rate quite high, with adaptive speciations 
happening every 1000 timesteps or so. If the mutation 
rate is too high, natural selection has no chance to weed 
out non-adaptive species, if too low, too much comput- 
ing resource is need to obtain interesting dynamics. In 
practice, the mutation rate is set about 2 orders of mag- 
nitude less than the critical amount needed for adapta- 
tion. In terms of speciation rates, the migration rate 
time constant might correspond to something of the or- 
der of 10 4 years, instead of the 10 years or so one gets 
from considerations of doubling times. 

This code is released as EcoLab 3.5. Due to a design 
flaw, performance of this code scales poorly with diver- 
sity, unless the code is run in parallel with one cell per ex- 
ecution thread. For this experiment, the runs took place 
on a 2 x 2 spatial grid, on a four processor parallel com- 
puter supplied by the Australian Centre for Advanced 
Computing and Communications, apart from one run of 
a 3 x 3 grid on a 9 processor system. Work is currently 
underway to implement a spatial version of the Ec<tab 
4.x code, which does not suffer from this performance 
problem. 

The results of a typical run is shown in hgure 6. The 
run started with a maximum migration rate of 0.01 at 
the bottom right hand corner of the figure and swept 
down to 10~ 10 before increasing. The migration rate 
was swept back and forwards 5 times over the 18 million 
time steps in the run. 

The first thing to note was that the expected response 
of diversity to migration rate was not there. We would 
expect a response of the form D oc A c , with A = 4 in 
the 2x2 case, and c varying smoothly between 1 for the 
infinite migration (panmictic) case and 2 for zero mi- 
gration. These results tentatively indicate that possibly 
c does not vary smoothly at all, but is nearly constant 
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Figure 7: Diversity growth on a 2 x 2 grid, with \~f\oo = 
1 x 10~ 6 . 

for most values of 7|oo. This needs to be resolved with 
further study. 

The second thing to note is the completely unexpected 
"resonance" at about 1 x 10~ 5 . It is not peculiar sta- 
tistical aberration, since the same result was obtained 
with completely different random number seeds, and fix- 
ing the migration rate at the resonance value produces 
an exponential growth in diversity (Figure 7). 

3 more tests were performed to determine if this re- 
sult is an artifact of discretisation, or a feature of the 
dynamics. The first involved changing the grid to a 3 x 3 
grid, which did not affect the location of the resonance. 
The second involved scaling all parameters in the model 
(r, /3, (j,) by 0.1, which is equivalent to changing the 
timescale. If the effect was purely due to dynamics, one 
would expect the resonance to shift one order of magni- 
tude higher on the scale, however little qualitative dif- 
ferent was observed. The third test involved performing 
the migration operator every 1000 timcsteps, instead of 
100. This did change the resonance value by 1 order of 
magnitude, ruling out certain classes of software faults. 

Conclusion 

The choice of diversity as a proxy measure for ecosys- 
tem complexity is a good choice. Complexity is obvi- 
ously constrained by diversity, so that bounded diver- 
sity dynamics also implies bounded complexity dynam- 
ics. However, in the case of evolutionary Lotka-Voltcrra 
dynamics, the system will tend to self-organise to a crit- 
ical surface where speciation is balanced by extinction. 
This surface defines the maximum allowed complexity 
for a given diversity value, which turns out to be pro- 
portional to the diversity. The analysis presented in this 
paper could be extended to other evolutionary ecologies 
as well. 

Whilst there is still debate about whether the bio- 
sphere is exhibiting unbounded complexity growth, I am 



persuaded by Benton's(2001) argument that the growth 
is nothing short of spectacular. In this paper I examined 
two possible mechanisms for diversity growth — speciali- 
sation which proves capable of delivering unbounded cre- 
ative evolution in ^ c tab, and biogeographic provincial- 
ism. Whilst I was only expecting biogeographic changes 
to deliver a modest impact on diversity, ^ c< tab deliv- 
ered a unexpected result of a "resonance" , where if the 
migration rate was tuned to this value, diversity grew 
exponentially. 
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